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Abstract. One of the main goals of design theory is to classify, characterize and count 
various combinatorial objects with some prescribed properties. In most cases, however, 
one quickly encounters a combinatorial explosion and even if the complete enumeration of 
the objects is possible, there is no apparent way how to study them in details, store them 
efficiently, or generate a particular one rapidly. In this paper we propose a novel method 
to deal with these difficulties, and illustrate it by presenting the classification of quaternary 
complex Hadamard matrices up to order 8. The obtained matrices are members of only a 
handful of parametric families, and each inequivalent matrix, up to transposition, can be 
identified through its fingerprint. 
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1. Introduction 

A complex Hadamard matrix H of order n is an n x n matrix with unimodular entries 
satisfying HH* = nl where * is the conjugate transpose and / is the identity matrix of order 
n. In other words, any two distinct rows (or columns) of H are complex orthogonal. Complex 
Hadamard matrices have important applications in quantum optics, high-energy physics [T] , 
operator theory [TB] and in harmonic analysis [TT], ^U\. They also play a crucial role in 
quantum information theory, for construction of teleportation and dense coding schemes [21] , 
and they are strongly related to mutually unbiased bases (MUBs) [8]. In finite mathematics 
and in particular, in design and coding theory, the real Hadamard matrices and their discrete 
generalization, the n x n Butson-type matrices [3] , composed of q-th roots of unity appear 
frequently. These matrices are the main algebraic objects behind the Fourier and Walsh- 
Hadamard Transforms and perfect quaternary sequences. They can also be utilized to obtain 
high performance quaternary error-correcting codes and low correlation Z4-sequence families 
[13], and shall be denoted through this paper by BH{q,n). The notations BH{2,n) and 
BH{4:,n) correspond to the real and quaternary Hadamard matrices, respectively. 

One of the main features of these BH{q,n) matrices is the large symmetry they have, 
which allows one to introduce some free parameters into them obtaining a continuous family 
of complex Hadamard matrices. This idea was used in a series of papers [5], [6], [7], [H], 
[18], where new, previously unknown complex Hadamard matrices were constructed in this 
way. The parametrization of BH{q,n) matrices is in some sense the continuous analogue of 
the switching operation, a well-known technique in design theory, and allows one to escape 
equivalence classes and therefore collect many inequivalent matrices into a single parametric 
family. 



Date: May, 2010., Preprint. 

This work was supported by the Hungarian National Research Fund OTKA-K77748. 

1 



2 FERENC SZOLLOSI 

The aim of this paper is to illustrate how powerful these continuous families are by pre- 
senting the complete characterization of the BH{4,n) matrices up to n = 8, and describing 
them in a compact and elegant form in terms of parametric families of complex Hadamard 
matrices. Interestingly, up to transposition, all equivalence classes can be identified through 
a recently introduced invariant, the fingerprint [17]. Constructing these matrices and iden- 
tifying a member from each equivalence class took only a few days of computer search, and 
we do not consider it a deep mathematical achievement. The main result of this paper, 
however, is the way we present the representatives of the equivalence classes, as with the aid 
of only a handful of parametric families of complex Hadamard matrices we can encode a lot 
of inequivalent BH{4:,n) matrices simultaneously. In principle, by introducing m free affine 
parameters into a BH{q,n) matrix one can obtain as many as inequivalent matrices of 
the same type, therefore describing an exponential number of equivalence classes through a 
single matrix. We believe that the concept of parametrization shall turn out to be a useful 
tool in discrete mathematics as well, and by utilizing this method one shall be able to obtain 
new lower bounds on the number of various combinatorial objects in the future. 

In the following section we give a full characterization of quaternary complex Hadamard 
matrices of orders 1-8, and describe them as members of parametric families of complex 
Hadamard matrices. Throughout the paper we adopt the notation from [TS] for well-known 
families of Hadamard matrices, such as F^^\a, b, c, d, e), etc. 

2. List of quaternary complex Hadamard matrices of small orders 

It is easily seen that quaternary complex Hadamard matrices of order n can exist only if 
n = 1 or n is even, and it is conjectured that this condition is sufficient [13]. It is trivial 
to check that there is a unique complex Hadamard matrix of orders n = 1,2, however, 
for composite orders starting from n = 4 one encounters infinite families already. It has 
been shown that all complex Hadamard matrices of order 4 are members of the following 
1-parameter family, depending on the unimodular complex number a given in [1] as 
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Note that all 2 x 2 blocks of H/^{a) are complex Hadamard matrices. Matrices arising in 
this fashion are called Difa-type (cf. [Hj). Recall that Hadamard matrices H and K are 
equivalent, if there exists permutation matrices Pi, P2 and unitary diagonal matrices Di, D2 
such that P1D1HD2P2 = K holds. We denote this property hj H ^ K. Clearly, one can 
consider matrices containing a full row and column of numbers 1 due to equivalence. Such 
matrices are called dephased. Let us remark here that the transpose of a Hadamard matrix 
is also Hadamard, but not necessarily equivalent to H . 
It is not apparent at first, how many equivalence classes of BH{A, 4) matrices are contained 
in the family H^^a) for various values of a. To quickly check that matrices -^4(1) and -f^4(i) 
are inequivalent, we can utilize Haagerup's invariant A{H), defined for complex Hadamard 
matrices H of order n, which is the following set [9]: 

A{H) := {hijhkihiihkj : i, j, k,l = 1,2, ... ,n] . 
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We readily see that A{H4{1)) and A(if4(i)) are different, therefore these matrices are in- 
equivalent. Also, by exchanging rows 2,4 one can see that H{a) ~ H{—a) and hence 
H{1) ~ H{—1) and -f^(i) ~ H{—i). Clearly the Haagerup invariant corresponding to 
BH{q, n) matrices can contain the g-th roots of unity only. 

Remark 2.1. It might happen that a complex Hadamard matrix H does not look like a 
BH{q,n) matrix, but is equivalent to one; however, BH{q,n) matrices can be quickly rec- 
ognized. Indeed, after dephasing H all of its entries should be some g-th roots of unity, 
otherwise A (if) would contain an entry which is not a g-th root of unity showing that H 
cannot be equivalent to a BH{q, n) matrix at all. This implies that checking the number of 
inequivalent BH{q, n) matrices in dephased parametric families is a finite process. 

The above remark and the paragraph just before it easily allow to show that there are 
precisely 2 inequivalent _Bif(4,4) matrices: -^4(1) and if4(i), respectively. 

An interesting 1-parameter family of complex Hadamard matrices of order 6 was discovered 
by Zauner [22] (and independently, but slightly later by Di^a [5]) who used the 2-circulant 
form of the following matrix 
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as a seed matrix to exhibit an infinite family of MUB triplets in C^. A quick computer search 
amongst all possible BH{4, 6) matrices confirms that despite the 1-parameter freedom we 
have here, there is only a single equivalence class of BH{4, 6) matrices. One might wonder 
why the matrices Dq{1) and DqIi) are equivalent. This is not trivial to realize, as one need 
to enphase one of the matrices to obtain the other one. This particular family of complex 
Hadamard matrices has been utilized in harmonic analysis as well. Considering the scaled 
logarithm of the BH{8,6) matrix Dq{{1 + i)/V2), one finds a 6 x 6 integer matrix having 
rank 3 over Zg leading to the disproval of Fuglede's conjecture in dimension 3 [TT]. 

Now we turn to the classification of BH{A, 8) matrices. An exhaustive computer search 
revealed after a few days the following 

Proposition 2.2. There are precisely 15 BH{A,8) matrices, up to equivalence; 5 of these 
matrices are equivalent to a symmetric one, while another 10 are inequivalent to their trans- 
pose. 

Now, having all inequivalent BH{4:, 8) matrices at our disposal, our goal is to describe 
them as members of certain parametric families of complex Hadamard matrices of order 
8. It turns out that all matrices found in Proposition 12.21 appear as members of 3 distinct 
well-known parametric families of complex Hadamard matrices that we will now present. 
For each of these families we will give the classes of quaternary Hadamard matrices among 
them. 
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The first infinite family we have here is the 5-parameter generalized Fourier family, stem- 
ming from the 8x8 Fourier matrix as the following: 



^8 {a,b,c,d,e) = 
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We have chosen the nonstandard parametrization here to emphasize the fact that one can 
obtain BH{4, 8) matrices stemming from the matrix above if and only if the values of a, b, c, d 
and e are fourth roots of unity (see Remark 12 .11) . 

Remark 2.3. As Fg^\a,b,c,d,e) features two rows and columns with entries ±1 it is an 
infinite family of jacket matrices (cf. [T2j). Note that Fg^\aw , bw'^ , cw^ , dw"^ , ew^) where 
w = is an eighth root of unity coincides with the matrix listed in [19], and hence is of 
Dita-type. The matrix corresponding to Fg^^(— i, — i/t^, — i, — i/t^) is equivalent to an 
infinite family of circulant matrices containing in particular Horadam's matrix K^lij (cf. [131 
p. 88]). 

After evaluating the matrix F^^\a,b,c,d,e) at the 4^ = 1024 possible quintuples, we 
obtained the following 

Proposition 2.4. The family Fg^\a,b,c,d,e) contains 8 inequivalent BH{A,8) matrices: 
the symmetric matrices F^^\l, 1, 1, 1, 1), Fg^\l, i, i, i, i), F^^\i, 1, i, 1, i), -F8^''(l, 1, i, 1, i) and 
further the matrices Fg^\l, 1, 1, 1, i), Fg^\l, 1, i, i, i), and their transpose. 

To prove that the matrices in Proposition 12.41 are indeed inequivalent, we need to recall 
a more powerful invariant than the Haagerup set. The fingerprint of a complex Hadamard 
matrix H of order n > A was introduced in [T7j and reads 

^H) := {{{v,{d),m,{d)) : t G I{d)} : = 2, . . . , [n/2j} , 

where I{d) is a finite index set counting the number of different moduli Vi{d) taken by the 
d X d minors of H, the multiplicity of these values being denoted by mi{d), respectively. 
Note that ^{H) is an ordered set. For instance, the fingerprint of the 8x8 real Hadamard 
matrix is given by ^{Fg^\l, 1, 1, 1, 1)) = 

{{(0, 336), (2, 448)} , {(0, 1344), (4, 1792)} , {(0, 1428), (8, 3136), (16, 336)}} . 

That is, for d = 2, 336 minors are and 448 have modulus 2, etc. In theory, the fingerprint 
is a more sophisticated invariant than the Haagerup set, because instead of looking at the 
2x2 minors only it takes into account the higher order minors as well, therefore capturing 
some of the global properties of complex Hadamard matrices. One of the limitations of the 
fingerprint is, however, that it cannot distinguish a matrix from its transpose. 

Now with the aid of this invariant one readily verifies that the 15 matrices appearing in 
Proposition 12.21 correspond to 10 different fingerprints, that is each matrix can be identi- 
fied, up to transposition, with its own. It is straightforward to check that the 8 matrices 
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appearing in Proposition 12.41 have 6 different fingerprints, and as the family F^\a,h,c,d,e) 
is self-cognate (cf. [19j) it follows that the transpose of the non-symmetric matrices are 

members of the family as well. Indeed, we have ^Fg^'*(l, 1, 1, 1, i)j = F^\l, 1, i, 1, 1) and 

(Ff (1,1, 1,1,1))"" = Ff(l,i,i, l,i). 

The list of the BH{4, 8) matrices is not complete yet, as in [H] another family of 8 x 8 
complex Hadamard matrices were obtained from tiling abelian groups: 
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It was shown in [H] that the matrix (1, 1, 1, 1) is not of Di^a-type, and hence a small neigh- 
borhood around it avoids the family Flf'\a,b,c,d,e) completely. The matrix corresponding 
to 5'8^''(1, 1, 1, 1) is equivalent to the "jacket conference matrix" Jg [I2]. By evaluating the 
matrix Sg^^ (a, b, c, d) at the fourth roots of unity we found the following 

Proposition 2.5. The family 5*8^^ (a, 6, c, d) and its transpose contain 8 inequivalent BH{4, 8) 
matrices: the real Hadamard matrix S'8^^(l, 1, 1, 1) and the matrix S'8^^(l, 1, 1, 1), which are 
equivalent to a symmetric matrix, and further the matrices 5*8^^(1,1,1,1), 5*8^^(1,1,1,1), 
5*8^'' (1, 1, 1, 1), and their transpose. 

By analyzing the fingerprint of the obtained matrices, one can see that the following 
equivalences hold: 5^^^(1, 1, 1, 1) ~ F8^^^(l, 1, 1, 1, 1) and (^5^^^(1,1,1,1)) -^8^^^(1,1,1,1,1), 

and therefore a further family is required to describe all of the 15 BH{A, 8) matrices. 
The following matrix was constructed from MUBs of order 4 by Dita [7J very recently: 
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We will see shortly that the family D^^{a,b,c,d,e) is essentially different from the families 
F^^\a,b,c,d,e) and Sg^\a,b,c,d), as it accounts for most of the BH{4,8) matrices. We 

remark here that the family (^D^^{a, b, c, 1, c)) is equivalent to the 3-parameter matrix Pg 
reported in [2]. We have the following 
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Proposition 2.6. The family Df^{a, 6, c, d, e) together with its transpose contain 11 inequiv- 
alent BH{A, 8) matrices, namely the real Hadamard matrix Dg^(l, 1, 1, 1, 1), Df^{l, 1, i, 1, 1) 
and the matrix Df^{l, i, i, 1, i), which are equivalent to symmetric matrices, and the matrices 
Dfsi^, 1, 1, 1, i), ^85(1' 1' 1' i' i)^ ^8b(1' 1' i' i' 1)' ^8b(1' i' i' i' 0; a'^t^ ^^e2r transpose. 

In particular, the matrix Dg^(l, 1, i, i, 1) is not a member of any of the families Fg^\a, b, c, d, e), 

iSg^^ (a, b, c, (i) or ^S'g^'' (a, 6, c, d) j . We are ready to state the main result concerning BH{A, 8) 
matrices. 

Theorem 2.7. There are precisely 15 inequivalent BH{4, 8) matrices, which can be described 
through the parametric families Fg^^(a, b, c, d, e), Sg^\a, b, c, d) and Dg^(a, 6, c, d, e). Five of 
these matrices are equivalent to a symmetric matrix, namely -^§^^(1, 1, 1, 1, 1), F^\l, i, i, i, i), 
Fg^^(i, 1, i, 1, i), Fg^^(l, 1, i, 1, i), and >S'g^''(l, 1, i, i). The other 10 matrices are given by Fg^^(l, 1, 1, 
Fg^^(l, 1, i, i, i), S'g'^^(l, 1, 1, i), 5'g*^''(l, i, 1, i), i5g^(l, 1, i, i, 1), and their transpose. These ma- 
trices, up to transposition, can be distinguished by their fingerprint. 

Finally, we would like to point out that another well-known invariant of EH {A, n) matrices, 
the Smith normal form (SNF, cf. [15J) can detect 6 equivalence classes only. In particular, 

the inequivalent matrices Fg^''(l, i, i, i, i) and iSg^^l, 1, i, i) share a common SNF. This is not 
surprising, as this invariant does not take into account the distribution of the minors of 
complex Hadamard matrices. 

3. CONCLUDING REMARKS 

In this paper we have described all quaternary complex Hadamard matrices up to order 
8. The number of equivalence classes for n = 1,2,4,6,8 read 1,1,2,1,15, respectively. 
By utilizing continuous parametric families we could present the inequivalent matrices in 
a compact and elegant form, which allows one to study, display, and store these matrices 
efficiently. We remark here that the classification of the nxn generalized Hadamard matrices 
over groups of order 4 has been completed up to n = 16 very recently [10], and we intend 
to continue this work and attempt to describe in a similar fashion the remaining BH[4:, n) 
matrices of size n = 10, 12 or even higher orders too in a forthcoming paper. 
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